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Abstract 

We comment on the recent paper by A. D. Dolgov, S. H. Hansen, S. Pastor 
and D. V. Semikoz (DHPS) [Astropart. Phys. 14, 79 (2000)] on the generation 
of neutrino asymmetries from active-sterile neutrino oscillations. We demon- 
strate that the approximate asymmetry evolution equation obtained therein 
is an expansion, up to a minor discrepancy, of the well-established static 
approximation equation, valid only when the supposedly new higher order 
correction term is small. In the regime where this so-called "back-reaction" 
term is large and artificially terminates the asymmetry growth, their evolution 
equation ceases to be a faithful approximation to the Quantum Kinetic Equa- 
tions (QKEs) simply because pure Mikheyev-Smirnov-Wolfenstein (MSW) 
transitions have been neglected. At low temperatures the MSW effect is the 
dominant asymmetry amplifier. Neither the static nor the DHPS approach 
contains this important physics. Therefore we conclude that the DHPS results 
have sufficient veracity at the onset of explosive asymmetry generation, but 
are invalid in the ensuing low temperature epoch where MSW conversions are 
able to enhance the asymmetry to values of order 0.2 — 0.37. DHPS do claim 
to find a significant final asymmetry for very large (5m^ values. However, for 
this regime the effective potential they employed is not valid. 
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I. INTRODUCTION AND OVERVIEW 



Active-sterile neutrino oscillations in the early universe can generate large differences in 
the number densities of relic neutrinos and antineutrinos .Q The resulting "lepton asym- 
metries" have important phenomenological ramifications: the suppression of sterile neutrino 
production because of large matter effects [|l|,@], a modification to primordial Helium syn- 
thesis if a sufficiently large asymmetry is created for the electron neutrino prior to or during 
the Big Bang Nucleosynthesis epoch 0], with also the possibility to describe inhomogeneous 
BBN scenarios IQ. 

We focus in this Comment on two flavour active-sterile oscillations. Furthermore we 
restrict our discussion to that region of oscillation parameter space which induces large 
lepton asymmetry creation well before neutrino decoupling at T ~ 1 MeV. An important 
issue is the "final value of the asymmetry", that is, the steady state value attained at the 
end of the dynamical evolution. The asymmetry is conveniently defined to be 

T — ^'^'^ ~ \ 
— 1 

where is the number density of species -0, and a is one of e, /i or r. In Ref. it is argued 
that the final value for L^^ is of order 0.2 — 0.37 for typical parameter choices. This claim 
is recently disputed by A. D. Dolgov, S. H. Hansen, S. Pastor and D. V. Semikoz (DHPS) 
in Ref. 0, who conclude that although the asymmetry can rise by five orders of magnitude 
above the baryon asymmetry level of ~ 10^^°, it cannot reach the 0.2 — 0.37 range found in 
Ref. [§. The purpose of this Comment is to explain the flaws in the DHPS analysis. 

The basic reason turns out to be very simple: DHPS neglect the pure Mikheyev-Smirnov- 
Wolfenstein (MSW) effect 0. As explained in Ref. [§, pure MSW transitions dominate 
lepton asymmetry evolution at lower temperatures, after the period of initial explosive Lj^^ 
growth. The approximate evolution equations derived by DHPS have sufficient validity to 
describe the onset of asymmetry growth, but they are not a good approximation after the 
explosive amplification phase ends. Their oversight of the role played by MSW transitions 
also leads them to exaggerate the importance of a "back-reaction term" [Bi in Eq. (51) of 
DHPS]. As we will explain, when this term is large, the MSW effect is also dominant so that 
their Eq. (51) is no longer valid. The size of the back-reaction is thus a moot point. 

We now clarify the structure of our paper, while simultaneously providing a summary of 
the main points to be made. These remarks will be of interest to those readers who want to 
know the gist of our argument but not the technical details. 

• Section |I| contains a quick survey of the Quantum Kinetic Equations (QKEs) govern- 
ing active-sterile neutrino evolution |§-TD]. These "almost exact" equations are not at 



issue; the debate is on their correct solutions. The QKEs simultaneously incorporate 



^For the historical record, the reader should note that Ref. |5| contains the claim that the neutrino 
asymmetry is always small (less than 10^'') and hence unimportant for all parameters of interest. 
This contrasts starkly with the later studies |l|-|3|. 
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coherent matter-affected evolution, decoherence due to non-forward neutrino scatter- 
ing, and repopulation of active neutrino distribution functions from the background 
plasma. 

Section |T| reviews how the essential physics encoded in the QKEs can be revealed. A 
fundamental issue is collision-dominated versus coherent matter-affected oscillation- 
driven amplification. The former pertains to conditions prior to and during the onset 
of asymmetry growth, in which case the evolution is well described by the static ap- 
proximation of Refs. [mH, and by the closely related adiabatic limit approximation 
developed in Refs. [|lT]-|T3[. Equation (51) of DHPS is identical to the leading order 



term of the static approximation if the Bi term is neglected, and is thus able to cor- 
rectly predict the existence of a critical temperature at which asymmetry evolution 
enters a brief explosive growth phase. Subsequent evolution, however, ceases to be 
collision dominated. In addition, a significant asymmetry, typically ~ 10~^ when ex- 
plosive growth stops, now exists in the plasma; the antisymmetric Wolfenstein term 
in the matter potential is sufficiently large to modify the neutrino and antineutrino 
oscillation patterns in noticeably dissimilar ways. A vitally important consequence 
of this is the separation of the neutrino and antineutrino MSW resonance momenta 
[0]. Depending on the sign of L,y^, one resonance momentum remains in the body 
of the (almost) Fermi-Dirac distribution, while the other is rapidly relegated to the 
tail. Figures lb and 2b depict this phenomenon for two different oscillation parameter 
choices. The separation means that the MSW effect, now dominant over the collision 
mechanism, has very asymmetric consequences for neutrinos and antineutrinos 0. For 
Lj,^ > 0, the resonance in the body is able to efficiently convert T^aS into z7s's, while 
the other is impotent in the scarcely populated tail. This important physics was not 
understood by DHPS. Indeed, when collisions are absent, neutrino density matrix evo- 
lution is driven only by matter-affected propagation (and the expansion of the universe 
of course). Equation (51) of DHPS does not reduce to the coherent matter-affected 
evolution equations, so it does not contain the pure MSW effect [|]. The static approx- 
imation, which yields ^^^^ = for zero collision rate, also neglects the MSW effect. 
Thus these equations are valid only in the collision dominated regime. Figures la and 
2a illustrate our arguments. The solid and dash-dotted lines represent respectively nu- 
merical solutions to the exact QKEs0 and the static approximation equations. Clearly, 
the latter describes the explosive growth phase excellently, but underestimates the 
subsequent growth of L^^. At steady state, the discrepancy is significant. The dashed 
curve reveals the physics behind the QKEs in a very amusing way. The generating 
code compares the decoherence or interaction length £int with the matter-affected os- 
cillation length im at every time step, with both lengths evaluated at resonance. If 
r|res = (Ant/^m)|res IS less than ouc, the static approximation equations are used to 
forward the system in time. For rlres > 1, pure adiabatic MSW evolution is employed. 
The quantity rlres therefore measures the relative importance of collisions and MSW 



^DHPS are unable to successfully integrate the QKEs, and do not accept these curves as accurate, 



although the numerical procedures used, and their stability, are thoroughly discussed in Ref. [14|. 
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transitions. In the event rlres -C 1, the MSW effect is damped since colhsions seriously 
disrupt the coherence. The opposite case has r|res 3> 1, implying that collisions are 
relatively negligible. The dashed curve in Fig. la agrees with the QKEs in both the 
high and low T regimes extremely well. There is a slight discrepancy immediately 
after explosive growth, since neither effect dominates at this stage — the very visible 
bump in Fig. 2a is an artifact of this "either/or" approximation in this intermediate 
regime. 

The DHPS equations underestimate the asymmetry after explosive growth even more 
seriously than does the static approximation, because of their Bi back-reaction term. 
Section |I^ compares Eq. (51) of DHPS with the static approximation. As already 
noted, and acknowledged by DHPS, these equations agree to leading order. We show, 
however, that the Bi term is large only when the MSW effect has taken over as 
the dominant driver of asymmetry evolution. Since DHPS Eq. (51) and the static 
approximation both neglect this effect, they should not be used when MSW transitions 
are important. (We will also identify where the MSW effect is unwittingly removed in 
their derivation.) Furthermore, the expansion parameter adopted by DHPS in their 
approximation scheme is not in fact the correct choice at lower temperatures when 
there is a significant asymmetry in the plasma. This misunderstanding led DHPS to 
erroneously claim validity for their equations at later stages in the evolution, especially 
the role of the back-reaction term which is but another casualty following from a faulty 
expansion parameter. In fact, the dramatic termination of the asymmetry growth due 
to this bogus term has its origin in the expression {l + x)~^ ^ 1 — a;. The left hand side 
corresponds to the static approximation (with a minor discrepancy), while the right 
hand side denotes the DHPS result, with x identified as the Bi term. The artificially 
strong cut-off arises because the right hand side is used even when x ~ 1. 

Section contains comments on other less important errors and some misleading 



statements made by DHPS, while Sec. ^ is a conclusion. 



The Appendix provides a translation of DHPS's notation into a more common lingua 
franca used by other authors. The rather obscure nomenclature of DHPS proves to be 
the biggest obstacle to understanding their analysis. 



II. SHORT REVIEW OF THE QUANTUM KINETIC EQUATIONS 

We consider the case of active-sterile two state mixing where the weak eigenstates z/^ 
{a = e, iJ, or r) and are linear combinations of two mass eigenstates z/^ and Uf,, 

Va = cos dQfa + siu Vs = - siu OoVa + COS (2) 

with 9q as the vacuum mixing angle, and cos2^o > by definition. The compact notation 

s = sin 26*0, c = cos 26^0, (3) 
and the convention Sm^ = ml — ml will be used from now on. 
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The one-body reduced density matrix [HHT^ for a z/q, ^ z/^ system of momentum p in 



the early universe can be parameterised through 

piy) = l[Poiy) + Piy)-a], (4) 

in which P{y) = Px{y)^ + Py{y)y + Pz{y)i is the "polarisation vector", and cr = a^^ + ayf + 
cr^z, where cxj are the Pauli matrices. The quantity ?/ is a dimensionless momentum variable 
that does not red-shift by definition: 

?/ = |ocp7^(^), (5) 

where TZ is the Friedmann-Robert son- Walker scale factor and T is the temperature. It is 
understood that the density matrices and the variables Pi{y) depend also on time t or, equiv- 
alently, temperature T, which, for me ^ T < m^, are related through dt/dT = —1/{TCT) ~ 
-Mp/5.44T^ where H = 7^/7^ is the Hubble parameter and Mp ~ 1.22 x lO^^ MeV is the 
Planck mass. 

We normalise the density matrices to give momentum distribution functions 

fuM = liPoiy) + PMKiy), UXy) = ^[^(y) - Pz{y)K{y), (6) 

for Ua and Ug respectively. The reference distribution feq{y) is of the Fermi-Dirac form [see 
Eq. (0)] with zero chemical potential. Similar expressions pertain to antineutrinos, with 
Pi{y) ~^ Pi{y)- The number densities are obtained from the distribution functions through 

= f^(y)y'^dy, (7) 

where ijj denotes the species. 

The evolution equations for P(2/) and -Po(2/) are HHH 



^ = V(;,) X P(;,) - D(y)[PM^ + P,(y)y] + i, 

The damping or decoherence function D is related to the total collision rate for of mo- 
mentum y with the background plasma, r(|/), via [^|TT],|T3[ 

D{y) = ^-^^\KGlT'y, (9) 
with ke — 1.27 and k^^ = k^- 0.92 (for rrie ^ T ^ ^^J.)■f\ It is useful to define the quantity 



^These values neglect the Pauli blocking factors. Including Pauli blocking factors lead corrections 
of order 10% for these quantities [p4|j6|. 
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-int 



= D-\ (10) 



to be identified as the decoherence or interaction length. 
The function feq{y) is the Fermi-Dirac distribution, 

where fla = fJ'Ua/T (and /i^ = jJ^v^T for antineutrinos) is a dimensionless chemical potential 
which does not red-shift but depends on the neutrino asymmetry. 

The matter potential vector P,pi?[, V(y) = Vx{y)^ + Vz{y)z, has coefficients 

v, = ?^. v; = -^^ + vi, (12) 

where io = 2yT/5m? is the vacuum oscillation length. The quantity Va{y) is the effective 
potential which takes the form 

V^{y) = ±v^Gpn^L(-) - ^Gpn^^^^, (13) 

where the +(— ) sign corresponds to neutrino (antineutrino) ensemble, Gp is the Fermi 
constant, Mw is the M^-boson mass, — 17 and ~ 4.9 (for mg < T < rn^), and 
L^"^ is the effective total lepton number (for the a- neutrino species). This last quantity is 
essentially a sum of the individual asymmetries: 

L(") ^ L,^ + L,^ + L.^ + L,^ + r/, (14) 

with 7] being a small term due to the asymmetry of the electrons and nucleons — typically 
of order ~ 5 x 10"^'^, calculated from 



7] = -^Ln (fora = e), t] = -^Ln (for a = /i,r). (15) 

For simplicity we shall henceforth write L instead of 

We now introduce Va and the dimensionless variables a and b by 

Va = ioVa = Ta + b, 



= -V2GFn^Lio ^ -8.033Ly (^j^^ 



VMeV 

6 /„^a2 



,.-V-2G,n,i,^ = -i-(r) fX_ . (16) 



where, in the definition of Va, —{+) pertains to neutrinos (antineutrinos), and 

T" ~ 18.9 (23.4) MeV, for a = e (/i, r), (17) 

is a convenient reference temperature. In this new notation, the matter-affected oscillation 
length is given by 
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^= = , . (18) 



We will also have cause to consider the ratio 

where 6 ~ 0.8 (2.0) x 10"^ for a = e (/i, r). In addition, a very important role is played by 



r = 



..nt VK^ + y.' _ v^^ + ,,,, 

D ~ ' ^^^^ 



which is the ratio of the interaction and matter-affected oscillation lengths. This quantity 
parameterises the relative importance of decohering collisions and coherent oscillations on 
the evolution of the density matrices. 

Note that the MSW resonance condition is given by = ^ c = f a = =Fo + ^, where 
— (+) is used for neutrinos (antineutrinos). Evaluating r at resonance, 



D 



(21) 



we find that the relative importance of collisions and coherent oscillations is strongly depen- 
dent on s = sin 26*0. 

To make our argument easier to follow, it is useful to distinguish two different asymptotic 
regimes for the resonance condition. 

1. Negligible total lepton number \a\ <^ \b\,c. This is the situation prior to the explosive 
growth phase, in which the resonance condition 

c = b, (22) 

is identical for neutrinos and antineutrinos. Equation ( p2D can be recast in the form 

rl(y)^r"(^)%-ici (23) 

where T^gg is the temperature at which neutrinos or antineutrinos of momentum y are at 
resonance. For the body of the Fermi-Dirac distribution, i.e., 0.1 < y < 10, one can 
then deduce from this expression the temperature range in which oscillations will be most 
significant. Note, however, that since Eq. (pTS]) holds only for T ^ 150 MeV, we are restricted 
to the regime |(5m^|/eV^ ^10^ for c ~ 1. 

Using the definition of T^^^ one may rewrite the function b via 

H=4^)'- (24) 

\ res / 

This will be useful later. Equation (|22|) can also be rewritten in the form 
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(25) 

which shows that the resonance momentum moves from low to high values as the temperature 
decreases. At y = y^gg, the condition \a\ <^ |6| = c is equivalent to 



1 1 



L « 1.0 (0.4) X 10-^ (^^j ' (y,",Jtct ^ L„ (26) 

for a = e(yU,r). The most relevant choice for y^^^ is about 2, because asymmetry growth 
starts approximately at this point as the resonance momentum moves form low to high 
values with decreasing T. This value is approximately correct if sterile neutrino production 
has been negligible until that moment (which is true for small enough mixing angles). 

If L ^ Lt, the Wolfenstein term a in the effective potential becomes dominant and 
two different resonances exist, one on either side of y^^^ in the momentum distribution. 
If L is positive, the resonance at the lower (higher) value occurs only for antineutrinos 
(neutrinos), and vice versa if L is negative. This "separation of the resonances" is crucial 
for understanding the mechanism of neutrino asymmetry generation. 

2. Negligible finite-temperature term \b\ -C \a\,c. This situation holds after a substantial L 
has been created. The resonance condition is now approximately given by 

c = Ta, (27) 

where — (+) refers to neutrinos (antineutrinos). In the case L > 0, the resonance occurs for 
antineutrinos (neutrinos) if Srn^ < (Sm"^ > 0). The opposite holds if L < 0. It is important 
to recall that the asymmetry generation mechanism actually requires < 0. 
Equation (^) implies that the resonance momentum takes the value 



tires Wres 



8.033 eV L(T/MeV)4 



(28) 



Notice that yl°^ < 7/°^^ when \L\ > Lt 



III. UNDERSTANDING THE LEPTON ASYMMETRY EVOLUTION CURVES 

The Quantum Kinetic Equations can be numerically integrated and lepton asymmetry 
evolution thereby extracted. In terms of the density matrices, the asymmetry is given by 



1 

2n^ 27r^ Jo 



Poiy) + Pziy) - Poiy) - Pz{y)\ fl{y)y'dy 



(29) 



Taking the time derivative of Eq. (|29|), and imposing a + s lepton charge conservation, i.e.. 



f 

Jo 



dPoiy) dPoiy) 



dt 



dt 



fea{y)y'^dy « 



dt \ 



0, 



(30) 
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one easily deduces that 



^ = — ^ / V,{Py - Py)flyHy. (31) 

In principle, one could use Eq. (129|) to directly calculate the asymmetry; Equation (^) is 
formally redundant. However, for very small asymmetries, the right hand side of the former 
equation is the difference of two numbers which are large relative to L^^. Equation (|3lD 
should therefore be used to circumvent the numerical stability problems potentially posed 
by the direct procedure. 

DHPS were apparently unable to achieve successful numerical integration of the QKEs 
for \5m?'\ > 10"^ eV'^. We would guess that at least part of the reason for this was that they 
did not use Eq. (^Tj). 

The solid lines in Figs, la and 2a show our numerically stable asymmetry evolution 
curves from integrating the QKEs. There are four phases in the evolution: Above the critical 
temperature T^, the asymmetry L^^ evolves from its small initial value to the t] ~ 10"^'' level 
such that L ^ 7]. At T = Tc, explosive growth occurs for a short period, during which the 
asymmetry rises to roughly 10~^ for the typical choice Sm^ ~ 1 eV^. After this phase, 
asymmetry growth continues at a more leisurely power law rate. One can measure from the 
graphs that L^^ oc T~^. Further evolution brings the asymmetry to a steady state value in 
the range 0.2 — 0.37 depending on the oscillation parameters chosen. 

The post explosive growth behaviour is disputed by DHPS. In this section, we review 
the physical explanation for the L^,^ curves obtained by solving the QKEs, and demonstrate 
concordance between these curves and those generated by certain approximate evolution 
equations which highlight the essential physical processes occurring in the various regimes. 



A. Collision dominated evolution and the static approximation 

The neutrino collision rate F increases as T^. We expect, therefore, that at high enough 
temperatures, collisions are sufficiently frequent to seriously disrupt the coherent matter- 
affected oscillatory evolution of the ensemble. The finite temperature b term in the effective 
potential also scales as T^, and predominates at high T if the initial neutrino asymmetry is 
small (say of the order of the baryon asymmetry). We focus first on this regime. 

A useful and intuitively reasonable measure of the importance of collisional decoherence 
is afforded by the quantity r\^^^ = (Ant/^m)lres defined in Eq. (^Tj). The collision dominated 
regime then corresponds to r\^^ <^ 1. In this regime, one can develop a physical picture 
for distribution function evolution that yields approximate time-development equations for 

and that reproduce the QKE behaviour exceedingly well. This technique has been 
called the "static approximation" IHQ]. From another perspective, the resulting equations 
have been shown to arise in the adiabatic limit approximation to the QKEs in the large 
damping or collision dominated epoch |ri|-pr^ . The approximate equations feature only the 



distribution functions, not the coherences Px,y, so we are working in a Boltzmann limit. 

We will not repeat the derivations in Refs. []I],0,|TT]-|r^ , but present just the results. 
Define: 
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ffa.s — fl^a, 



^ eq ^ eq 

The z's are the actual z/q, and Us distribution functions normahsed by the Fermi-Dirac 
distribution with zero chemical potential. In the static approximation, the sterile neutrino 
z^s obey the rate equations 

dz T ffz T 

'^^s ^ as r 1 '^^s as r -i ^QQ"\ 

'dT ~ ~T-cr ~ ' 'dT ~ ~l-cr ~ ' 

where 

_ 1 Ts^ _ 1 Ts' 

°' ~ 4s2 + d2 + (c-6 + a)2' - 4 s2 + cf2 + (c-6-a)2' ^ 

are the transition rates. 

Using lepton number conservation and the following definitions: 



-L ~t~ ^S — ^S ■^S 

zj = , Z„ 



2 ' " 2 ' 

we see that the neutrino asymmetry evolution equation is given by 

dL„ 1 foo dz 



(35) 



dt n^2TT^Jo dt ■"''^ 



' f'y'dy. (36) 



From Eq. (Rl) the rate of change of z^ is obtained by solving the coupled equations 



dz~^ F + F F — F 

dz~ F — F F + F 

'^'^s ^ as ^ as r 4--! , as J- c 



z:-zT + 



dt 2 ' " 2 

Adopting the instantaneous repopulation approximation 

1 + ey 1 + ey 



Z 



(37) 



(38) 



and expanding these expressions to first order in the chemical potential, we finally obtain 

dL _ f-o AL + BL,^ + B' 
d{T/MeV) ~ Jo A{s^) ^^^^ 

with 



A ~ -1.27s^A;, 



5m2 / VMeV 



) {b-c)y'il + eyr\l-zt), 



B ~ 0.09s2fc„ f-^V + d^ + ic~ hf + a^]y^ey{l + e^)-^, 
ViVleV / 

B' ^ -QMs^K (^) [s^ + d^ + {c- hf + a']y\l + e^)"^^;, (40) 

where we have used = 2^(3)T^/7r^ ~ T'^/4.1, and 

A{s^) = [s^ + d^ + {c-b + af] [s^ + d^ + {c-b- af]. (41) 

Equations (|37|) and (|39D form a coupled system that must be simultaneously solved to track 
asymmetry evolution and sterile neutrino production. 
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B. Repopulation, the damped MSW effect, and the adiabatic Umit in the coUision 

dominated epoch 



References [|lT|-p!^ show how the static approximation follows from the QKEs in the 
adiabatic limit. These results are briefly reviewed, with the focus on repopulation handling. 
We will see that repopulation is not a critical issue for asymmetry evolution in the collision 
dominated regime, although it does play an interesting role in the otherwise mundane T > 
epoch. Importantly, we explain how the static approximation excludes the MSW effect. 

One may extract the collision dominated adiabatic limit beginning with the four coupled 
QKEs [Eq. (^] [|T^. These calculations show that the proper incorporation of repopula- 



tion affects the evolution of L^^ at sub-leading order only; an excellent description of the 
asymmetry evolution in the said limit may still be obtained from less rigorous repopula- 
tion treatments. We will not repeat here the somewhat involved algebra of Ref. [0, but 



instead explain the gist of the result through a simpler argument. The key is to consider 
repopulation in two extreme limits: a full decoupling limit and a thermal equilibrium limit. 
The first limit corresponds to completely switching off the refilling interactions, i.e., ^ = 



|Tl|], while the other assumes f^^ = /eq, where the right hand side includes the appropriate 



chemical potential [Q. In comparison, the latter limit is much closer to the truth since 



weak interaction induced scattering occurs frequently in the epoch under consideration. 
In the full decoupling limit, the QKEs [Eq. (||)] reduce to the homogeneous system |jlT 



9P 

^^KP, (42) 

where the 3x3 matrix /C is given by 

-D -V, \ . I -d {c±a-h) 
fC= \ V, -D =- -(c±a-6) -d -s\. (43) 

/ ^0 V s 

Rewriting Eq. (14^) in the instantaneous diagonal basis, and taking the adiabatic limit by 



ignoring the time derivatives of the diagonalisation matrix |TT|-|T3[, the evolution is now 
driven by the eigenvalues Ai,2,3 of K. In most cases, Ai_2 form a complex conjugate pair. 
Since we are interested in the collision dominated regime which includes requiring that 
VxjD <^ 1, these are well approximated by 



Ai,2 i^V^ + = -irl ± (44) 

The remaining eigenvalue A3 is small, real and negative, and obeys [1^ 



Dv^ Dv: 



X „, ^ • X 



D-^-U-uT- (45) 



' {D-\,f + V^ + V^- D^ + V^ + V^ 1 + ^' 

Since the real parts of Ai,2 are negative and much larger in magnitude than A3, the solution 
lies in the direction of the eigenvector associated with A3, while the transverse components 
are strongly damped (which is another way to define collision dominance 0,0). 
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It is very important to observe that the oscillatory aspect of the evolution, driven by the 
imaginary parts of Ai^2; is severely damped in the regime concerned. This is precisely where 
the pure MSW effect is removed from the dynamics. The DHPS procedure also contains 
a step which nullifies MSW transitions, a point they failed to notice. Note that the ratio 
of the imaginary to the real parts of Ai,2 is approximately tint/C-m- When this ratio is not 
small, many oscillations happen within one interaction length, thus allowing coherent MSW 
transitions to take place, particularly at low temperatures where £int is large. Using the 
diagonalisation matrix, one can also see that the oscillation amplitude is large. See Refs. 
| jTT| , P^ for further discussions on extracting the MSW effect from the QKEs. 

The sterile neutrino production equation obtained in the fully decoupled limit is simply 

^--\Hzc.-Zsi (46) 

which is in agreement with the static approximation. The related asymmetry evolution 
equation is also of static approximation form. 

In the opposite thermal equilibrium limit, the QKEs become also a homogeneous sys- 
tem of three differential equations if chemical potentials are ignored in the active neutrino 
distribution functions, so that ^ = and thus ^ = — with the matrix K, given by 



/C= V; -D -VA=-\ -{c±a-h) -d -s\. (47) 





Repeating the same procedure as before yields only one difference: the denominator of 
A3 now requires the substitution V^^ — > This, however, has no practical importance 

in the collision dominated regime, since the condition r\^^^ = (^int/-^m)lres ~ (^/-^)lres ^ ^ 
is required to validate the approximate equations at resonance. As r\^^ — > 1, collision dom- 
inance collapses locally; MSW transitions become the chief amplification mode P,P|,pT}^TB[] . 

To conclude, we find that the rate of sterile neutrino production and asymmetry genera- 
tion is the same at leading order in V^, independent of the description of the active neutrino 
distributions: full decoupling thermal equilibrium |]T^, or the more general equation 



T3| . Deviations from the initial thermal equilibrium distribution have no noticeable ef- 
fects on the asymmetry generation mechanism. This statement is qualitatively understood, 
because strongly damped oscillations do not produce significant spectral distortions. 



C. Comparison of the QKEs and the static approximation 

The dash-dotted curves in Figs, la and 2a show the integration of the static approxima- 
tion equations (pTf ) and (pQ]), which are in excellent agreement with the QKE results (solid 
lines) prior to and during explosive growth. The essential physics in the collision dominated 
regime is very well captured by the static approximation: the MSW effect is damped, and 
rate equations for the neutrino distribution functions govern the evolution. In terms of the 
QKEs' adiabatic limit, the transition rates in these equations are simply the eigenvalue A3. 

Moreover, one can easily infer from Eq. (|39D the existence of a critical temperature Tc 
and its approximate value by examining the dominant term A P|J^JTT[] . This term has a fixed 
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point at L = that is stable above Tc, and unstable below it where deviations from L = 
will cause runaway positive feedback. The transition is controlled by the function {b — c) in 
Eq. (^Ol). Since b ~ T^, the A term is generally positive for < at high temperatures, 
and L is destroyed. At T^, the "effective A" [the y integration in Eq. (|39D must be performed 
to determine the net effect of A] becomes negative, leading to a spurt of exponential growth. 

The subdominant B and B' terms in Eq. ( ^Qf ) control the extent to which the condition 
L = is approximated prior to Tc] realistically L = does not exactly define a fixed point 
since B and B' are not proportional to L. Numerically, we find that L ~ 10~^^ is a typical 
value. Interestingly, DHPS report values as low as 10~^°°. The reason for this is clear: 
the B and B' terms in Eq. (|39| ) arise from the nonzero chemical potentials in and Zq. 
These terms are absent in the DHPS equations since chemical potentials are neglected in 
the repopulation term, even when they attempt to numerically solve the QKEs [see Eq. (9) 
of DHPS and following discussion]. In their case L = is much closer to being a true fixed 
point. This is a relatively minor issue, and not germane to the controversy over the final 
asymmetry value. Nevertheless, it is of some numerical interest since tracking values of order 
10~^^, rather than 10"^°'', is a much easier task. Also, if 10"^°'' was the correct number, one 
must consider the presence of statistical fluctuations in the early universe as they would be 
able to change randomly the sign of the asymmetry at different spatial points |P]J5|JT^. In 
the presence of chemical potentials in the active neutrino distributions, however, this effect 
can be excluded |]14[ for a large range of mixing parameters (roughly for sin^26'o ^ 10~^). 

On the practical front, the neglect of chemical potentials when numerically integrating 
the QKEs leads to an artificially higher sensitivity to small effects on the solution (numer- 
ical error, oscillatory terms). DHPS observed this sensitivity, but did not realise that the 
inclusion of chemical potentials was remedial from this perspective. This may be one rea- 
son why DHPS found that their results were chaotic in sign for some mixing parameters, 
while an analysis performed with chemical potentials did not |T^. DHPS took this synthetic 
chaoticity as grounds to not rely on their own numerical solutions to the QKEs, or those of 
others 1,0.0 

After explosive growth. Figs, la and 2a show that the static approximation equations 
significantly underestimate the subsequent augmentation of L. This is the epoch in which 
the adiabatic MSW effect takes over as the principal propellant. 



D. MSW dominated regime 

Collision dominance and thus the static approximation collapse when iint/^m approaches 
and becomes larger than one MM- This is particularly evident in the adiabatic limit approach 



11,12]. Consider the eigenvalues of the matrix K, at resonance: 



■^Curiously, DHPS relied on the results of Ref. [17], which also reported rapid oscillations in the 
sign of the neutrino asymmetry for dm"^ < eV^, although the equations used in Ref. 1 17] are no 
more accurate than those in Refs. []^Jl^. DHPS may not have been aware of this incongruity since 
Ref. |]l^] does not contain any figure of the asymmetry behaviour, although rapid sign changes are 
clearly reported. 
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ALs = -A • (48) 

When collision dominance breaks down, i.e., iVx/D)\j.^^ > 1, the oscillatory imaginary parts 
are not preferentially damped. The real eigenvalue A3 and the real parts of the complex 
conjugate pair are now comparable in size. The imaginary components drive the MSW 
effect [|TT],|T2[, inducing many large amplitude oscillation cycles within the decoherence time 
scale. 

The simplest way to explore the implications of the MSW effect is through the approach 
of Ref . . We begin with a very important observation: after the initial burst of L creation, 
the neutrino and antineutrino resonances become separated in momentum space. This is 
illustrated in Figs, lb and 2b, and is easy to see analytically from Eqs. (0) and (pTj). At 
Tc, Eq. (|22|) holds, with neutrinos and antineutrinos having the same resonance momentum. 
After a substantial L has been created, the resonance condition evolves to being given by 
Eq. (pTl) , and the resonance momenta are separated. If L > 0, then the neutrino resonance 
quickly moves to the tail of the distribution, while the antineutrino resonance stays in 
the body. The MSW effect therefore has very asymmetric consequences for neutrinos and 
antineutrinos: the former are unaffected, while the latter continue to be processed into sterile 
states, thereby adding to the asymmetry. We shall henceforth assume L > for definiteness. 

Let us focus now on the large region of parameter space for which the MSW transitions 
are adiabatic. Here, complete conversion of antineutrinos at the resonance momentum into 
sterile states means that the asymmetry is augmented, while the position of the resonance 
evolves with L. Indeed, the rate of asymmetry growth depends on how quickly the resonance 
moves through the distribution. In the case of a small resonance width, it is easy to make 
this connection heuristically, which in turn gives rise to a simple evolution equation |^ 

dLy^ IT , , , , W 2 ^Vres _ -jr^C^Z/rcs ( AOi\ 

Using the resonance condition of Eq. (pT]), one obtains from Eq. p9| ) the expression 

dL^^ 4 Xy^es 



dT T 1 + 4^ 



(50) 



as a non-linear asymmetry evolution equation, for the case where < 0. Repopulation is 
incorporated as per the procedure described in Ref. |jT8|. 



The result of solving Eq. (|50D was first presented in Fig. 1 of Ref. [Q, and compared with 
the QKE solution for temperatures T ^ Tc/2. Excellent agreement was obtained, and the 
final asymmetry was found to be in the 0.2 — 0.37 range depending on the oscillation param- 
eters. Numerical solutions of the QKEs and of Eq. (|5^ yielded completely consistent results 
for the temperature range where the latter equation was expected to be a good approximation. 

Notice in particular that when L^^ <^ 1, Eq. (|50|) simplifies to 

dL,. ALy , ^ 
^ (51) 



dT 



1 
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which immediately imphes that L^^ oc in that regime. The behaviour that can he 
empirically measured from the QKE curves is fully explained by the adiabatic MSW effect. 

As the asymmetry continues to grow, the antineutrino resonance evolves out of the 
body into the high momentum end of the distribution. The asymmetry becomes frozen at 
some final steady state value, whose approximate magnitude can be easily understood by 
integrating the Fermi-Dirac distribution from y = yl°J ~ to ?/ = oo: 

.final J_ (T^y rjM^^ 3 (T^y 

^ 4C(3) \ T^) Jo l + ev 8\tJ ' ^ ^ 

The temperature ratio takes care of reheating due to e+e~ annihilations at T ^ rrif. ~ 0.5 
MeV. Equation (|52|) gives the approximate magnitude only, because the distribution changes 
with time as the asymmetry is created. Numerically, the final values found by incorporating 
proper thermalisation effects are quite close to this estimate. 



E. Combining the collision dominated and MSW dominated regimes. 

Figures la and 2a summarise neatly the overall picture of asymmetry evolution, as dis- 
cussed in Sec. |I[ The solid and dash-dotted lines represent solutions to the exact QKEs and 
the static approximation equations respectively. Clearly, the latter offers an excellent de- 
scription for the explosive growth phase, but underrates the subsequent growth of Li,^. The 
dashed line is generated by a code that computes the ratio between the resonance interaction 
and matter-affected oscillation lengths, r|res = ( Ant /^m) Ires' every time step. If r|res < 1, 
the static approximation equations are used to advance the system. For rl^es > 1, pure 
adiabatic MSW evolution is employed via Eq. (|50D. In both the high and low T regimes, 
the dashed line and the QKEs are in superb agreement. The small disparity directly after 
explosive growth arises from the fact that neither effect is dominant in this intermediate 
stage. We conclude that the physics of asymmetry growth has been thoroughly understood. 



IV. COMPARISON WITH THE DHPS RESULTS 



In Sec. [IV A] we "deconstruct" the DHPS equations, while we devote Sec. [IV B| to pin- 
pointing where in their derivation the MSW effect is unknowingly neglected. 



A. An analysis of the DHPS equations 

For comparison with the static approximation, we rewrite DHPS's evolution equations 
in the notation of this Comment using the "conversion relations" given in the Appendix. 
Equation (42) of DHPS describes sterile neutrino production. Its translation is 



d{T/MeV) AHT A(0) 



Comparing this with Eq. ([37[), one sees that it agrees with the latter's leading term, when 
(i) chemical potentials are neglected and thus z+ ~ 1, and (ii) the replacement 
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A{s') A(0), 



(54) 



is made. 

The neutrino asymmetry evolution equation [Eq. (51) of DHPS] is equivalently 



1 



dL 



L d{T/MeV) Jo A(0) 



A 



1 - 



g2 d^^(c- by + 
4 A(0) 



dy. 



(55) 



The first term on the right hand side is again connected to the leading A term of the static 
approximation [Eq. (|39|) ], with the substitution Eq. (^If). The B and B' terms are absent as 
expected, since DHPS neglect chemical potentials in the z/^ and z/q, distribution functions. 
The second term in the square brackets is DHPS's back-reaction term Bi. This suppos- 



edly new term has, in fact, an interesting interpretation. 
One can perform the expansion 

d' + ic- bf + 



Consider the function A(s /8). 



A(s2/8) 




A(0) 



+ {c-b + af' + {c-b 



(56) 



Thus the back-reaction term is but the first order correction in an expansion of A{s^ j'S). But 
we have already discussed in Sec. [ill ij| this sort of expansion, in the context of the V!^ term 
in the denominator of the eigenvalue A3, where the full decoupling and the instantaneous 
repopulation limits yielded A{s'^) and A{s'^/2) respectively. We argued that this distinction 
was of no practical importance, because the numerical difference between these two cases 
was only significant when the collision dominance assumption used in their derivation broke 
down. In a nutshell: (i) The first order terms are only large when the MSW effect is 
important, (ii) When the MSW effect is important, equations such as obtained from the 
static approximation, and DHPS (51), are no longer a good approximation to the QKEs. 
Let us now evaluate Eq. (^) for the antineutrino (c — 6 — a = 0) resonance: 



(57) 



When the asymmetry grows at the critical temperature Tc, the resonance condition passes 
from a regime in which 6 ~ c when L is negligible, to a ~ c. In this second regime the term 
<^ and the last expression becomes simply. 



A 


A 






A(sV8) 


res ^(^) 


res 


[ 4rfLKes + 4a2) 



rf2 

""res 



A 


A 






A(sV8) 


res ^(^) 


res 


[ 8<,J 



(58) 



So we see that the expansion parameter on resonance is s'^/d'^^^, i.e., the square of the ratio 



between the local interaction and matter-affected oscillation lengths 



As explained 



earlier, a small r\j.^^ denotes that collisions are frequent enough to interrupt coherent evo- 
lution and thus to prevent the MSW effect, while a large r\^^^ means MSW conversion can 
occur as the coUisional mechanism switches off. Of course this transition is gradual. 
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Moreover, using the equality d = Sb, one can also see, using Eqs. ( p4D and (|25|), that 



res 




(59) 



res 



where y^^^ is the resonant momentum when the asymmetry is negligible. Recall that |/° 
is very small at high T, but grows as the temperature decreases, reaching ~ 2 at T^. The 
quantity i/res is the general resonance momentum, which attains a minimum value depending 
on Sm^ (smaller for low 5m^) when the asymmetry grows (see Figs, lb and 2b). It is thus 
clear from Eq. (|59D that if one continues to use Bi in the evolution equation when the 
asymmetry is large and /d"^^^ ~ ^^(1); the rise of the asymmetry will be artificially made 
slower. From that moment its evolution will be approximately \L\ ~ 50 Lt and from the 
equations Eq. (|26|) and Eq. (|25|) this implies that \L\ oc 1/T, a behaviour observed by DHPS 
in their results. In this way the solution can grow just one order of magnitude since the 
onset of the instability. 

In summary, the back-reaction term Bi is generated by an expansion in s/dj-es, and not 
simply s/6 as DHPS believe. Their evolution equation is therefore invalid for temperatures at 
which s^/d^es > 1- As emphasised before 0,0, both the static approximation and the DHPS 



equations belong to a family of collision dominated idealisations which require s/c? -C 1 (plus 
other conditions) for their validity. When s/d > 1, the MSW effect is missed in both cases. 
Clearly DHPS did not properly appreciate the constraints on the real expansion parameter 
adopted in their derivation, and hence trusted their equations' validity for all temperatures. 

For very large 5m^, DHPS do actually observe a large final value. This is because 
|L|fin ~ lQ-^{\5m'^\/eV^y/^{TjTf) ~ 10-^{\6m'^\/eV^y/^ where Tj is the freezing tempera- 
ture corresponding to a resonant momentum (of antineutrinos if L is positive, of neutrinos 
if L is negative) well in the tail of the distribution (y/^g > 10) and is very weakly dependent 
on the exact choice of y{^^. Thus for |5m^| ~ 10^ eV^ they find |L|fin ~ 0.1, and are there- 
fore still able to state that active-sterile neutrino oscillations can yield large asymmetries. 
However, this high brr? asymmetry rise is spurious. For dm? ^ 10^ eV^, the effective poten- 
tial as given by Eq. ( |7^ is meaningless, since the critical temperature is now much larger 
than 100 Me V. The conditions of the plasma are different, especially at temperatures above 
the quark-hadron phase transition. The dynamics of neutrino oscillations in this regime is 
actually an open problem. 



B. How DHPS neglect the MSW effect 

In extracting their simplified evolution equations from the QKEs, DHPS make use of 
two small expansion parameters: (i) 1/Q, and (ii) sin 2^. We now examine their merits. 

1. Small 1/Q expansion. The derivation begins with an attempt to solve analytically two 
out of a system of eight coupled differential equations, namely, Eqs. (27) and (28) in DHPS: 

h'_/Q = ui^ - 7/i_ - 

l'_/Q = L{^a_-s_)-Uh_--il_ + VZh+, (60) 
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where prime denotes differentiation with respect to r. These expressions are equivalent to 



DHPS devise an oscillation averaging procedure that consists of setting all terms containing 
1/Q (i.e., the r-derivatives) to zero, from which they obtain what are actually "steady state" 
solutions to the variables /i_ and We now explain why this method is flawed. 

Firstly, the variable Py enters the big picture as the "driving force" for Vs production 
^ [see Eqs. (|3T|) and (|36|)]. In the context of solving differential equations, one may ignore 
contributions from oscillatory components in the driving force Py only if the oscillation 
amplitude is much smaller than the size of its steady "main term". Otherwise, oscillations 
propel the evolution, producing a marked change in the integrated variable Zg (i.e., the 
MSW effect) when the oscillation frequency is a minimum at resonance, in the same way 



In their derivation, DHPS do not consider at all the magnitudes of these oscillatory terms, 
and discard them on the basis that they impair computational efficiency, when in fact even 
the adiabatic solution to Py (and Py) is purely oscillatory at resonance in the post explosive 
growth environment where a relatively large r|res = (^mt/^m)|res ratio prevails Their 
equations cannot describe the sharp change in Zs at the low temperature resonance because 
the dominant and strongly frequency- dependent term responsible for the MSW effect has been 
removed. Indeed, one cannot know the sizes of the steady and oscillatory components in Py 
or /_ without some input of the initial conditions; to eliminate dynamical variables in the 
equation's differential form is premature and may be dangerous.^ 

Moreover, the process of artificially factoring out a large quantity by redefining the 
integration variable and then exploiting its inverse as an expansion parameter is not in fact 
a well-justified manoeuvre to obtaining steady state solutions. To see this, one needs only to 
consider the case of exponential growth/decay, ^ = Ax, where the approximation — 
and thus the conclusion x ~ crumble for a positive A. 

Fortunately, steady state solutions to h_ and Z_ are indeed valid for some oscillation 
parameters {5m? , sin2^^) when collisions are frequent enough to quickly damp the oscillations, 
even at resonance, provided that the adiabatic condition is also satisfied |jTl],|l2|. Thus DHPS 
have managed to reproduce results (albeit from a somewhat dubious procedure) that are 
consistent with the static approximation, in the manner that the latter may also be derived 
in the d ^ s limit by setting ^ ~ ^ ~ 0. But, like the static approximation, their 
equations do not apply to the low temperature region after explosive growth where steady 
state solutions do not hold at resonance. 



^We also mention a point noted in Ref. that when the asymmetry has grown sufficiently to 
separate the neutrino and antineutrino resonances, the differences in the distributions, and their 
derivatives, can be as large as their sums, simply because one distribution is on resonance while 
the other is far from it. 





18 



2. Small sin 26 expansion. This approximation is used on two occasions, first as a pretext 
for disregarding tlie variable Pz — Pz (or a_ — s_) in Eq. (|60| ) [Eq. (28) in DHPS], and tlien 
in evaluating the eigenvalues and other properties of the matrix Ai in their Eq. (33). 

The first exploit is precarious from the perspective of solving differential equations, since 
Px,y,z,o (and therefore a^-, s±, h±,l±) are interdependent dynamical variables and one has 
no prior knowledge on their sizes. Indeed, it is far more seemly to assume a priori the 
Pz — Pz term to be of some substance at resonance when the MSW effect kicks in at low 
temperatures, given that only antineutrinos (for L > 0) can be converted into sterile states 
in this epoch by the said mechanism. The neglect of this term is also related to DHPS's 
observation of the asymmetry reaching arbitrarily small values prior to explosive growth, to 
be discussed later in Sec. 

In the second instance, the expansion procedure employed to diagonalise the matrix 



M 



/ F \ 

-27 -F 
-7 f/ 



U = U\1- 



72 + f/2 



7 = 71 + 



72 + f/2 



(62) 



where F = —s, 'j = S/r"^ = d, D = —a oc L and U = l/r^ — 1 = 1 — 6, is not itself at flaw. 
However, it does have its own limitations which DHPS did not notice. 

Equation ( |62D suggests that the term F (or s) should be compared with 7, U, 7 and U 
if it is to be treated as a perturbation. DHPS consider only s and S, presumably because 
7 = (5 at resonance, which they define as occurring at r = 1 such that U = (or c = b). 
However, this resonance condition assumes a small asymmetry, i.e., |6| ^ \a\. Thus it is able 
to describe the high temperature regime where the term b dominates, but not after a sizeable 
asymmetry has been created, where the now separated neutrino and antineutrino resonances 
occur at any two t's other than r = 1. The ratio s/6 is meaningless in this second regime. 

Lastly, it must be stressed that MSW flavour conversion is not a perturbative effect. 
Small sin 2^ expansion equals no MSW effect, in the spirit that treating the off-diagonal 
sin 26' in the Schrodinger Equation for a two-flavour coUisionless system as a perturbation 
will not yield the correct behaviour at resonance. 



V. OTHER COMMENTS 
A. The mystery of the oscillation in sign 

DHPS allege that results presented by other independent groups are in conflict on the 
issue of oscillations in the sign of the asymmetry at the critical temperature. They use this 
to cast doubt on the credibility of some of these works. The true status of this subject is 
actually very clear, as we will explain below. 

To begin with, one must realise that two different systems are discussed in the literature. 
The flrst is a simplifled scenario in which all neutrinos have the same momentum, usually 
taken to be the zero chemical potential thermal average, p = (p) ~ 3.15T. The second is the 
realistic case where the neutrinos are described by a Fermi-Dirac distribution. The mean 
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momentum simplification should be viewed as a "toy model" which is not always a good 
approximation to the realistic case. However, it does give some useful qualitative insight. 

Another critical issue is whether the exact QKEs or an idealisation such as the static 
approximation are used in the calculations. 

Using the mean momentum approximation, Ref. observed oscillations of sign for 
some parameter points and apparently drew the incorrect conclusion that oscillations were 
a generic feature of lepton number generation. The simplified system was revisited in Ref. 



2T| , where oscillations were found to occur for only part of the parameter space (see Fig. 2 
of Ref. 1^). In both papers, the mean momentum QKEs were solvedQ 

In the realistic Fermi-Dirac case, the issue of the sign of L,^^ was first discussed in Ref. 
in the context of the static approximation; in the region of parameter space where this 
approximation is valid, the sign does not oscillate. However, as also pointed out in that 
paper, the static approximation breaks down at large sin^ 29 due to the rapid rise in lepton 
number during the explosive growth phase. Similar results were obtained in Ref. ||23|| . 

The oscillation issue was addressed in detail in Ref. , where comparisons of the static 
approximation and the QKEs showed that they gave concordant results for sin^ 29 ^ 10~^. 
For larger values, sign oscillations were found using the exact QKEs, but not the static 
approximation simply because the latter is invalid for too large a sin^ 29. Also, Refs. |n 



and [£1[] are not in mutual disagreement, since the former examined the realistic case with a 
Fermi-Dirac distribution, while the latter considered the mean momentum approximation. 



B. Abrupt cut-off at large mixing angles? 

In an attempt to locate the region of parameter space for which neutrino asymmetry 
generation is possible, DHPS found the following approximate upper and lower bounds on 
the so-called "region of instability" : 

sin^ 29o ■ ^\6m^\/eV^ < 10-^ (63) 
sin2 2eo ■ dWl/eV')^ > 4 x 10"^^ (64) 

While the second constraint agrees with that found in Ref. the first is too severe compared 
with previous calculations P,p3| .p] In addition, DHPS observe an abrupt transition from the 
"unstable" to the "stable" (i.e., no growth) regions. This is also at odds with earlier works, 
and is most likely a byproduct of numerical errors for the following reason. 

Increasing sin 26' for a fixed Sm"^ intensifies sterile neutrino production, which, in turn, 
has the effect of shifting the critical momentum to values larger than 2 |0. By Eq. (p3[). 



^ Recently a detailed analytic study of the mean momentum model was presented in Ref. . 

^The lower bound in Ref. Q is superficially more stringent than Eq. (^) (5 x lO"^'^ instead of 
4 X 10~^^) simply because a more restrictive condition — that the asymmetry must grow to the 
~ 10~^ level — was imposed in Ref. |^], while DHPS demand only amplification by as little as one 
order of magnitude. However, both generate the same slope on a (log (^m^, log sin 2^) graph. 
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the critical temperature decreases correspondingly, and explosive asymmetry generation is 
thereby delayed. (The growth curve also becomes smoother but a substantial asymmetry is 
still eventually attained.) It follows that for very large mixing angles, no explosive growth 
can occur since the critical momentum is then banished to the tail of the distribution. 

DHPS note correctly the connection between suppression of asymmetry amplification for 
a large sin 29 and a substantial sterile neutrino distribution manufactured prior to the onset 
of an otherwise explosive growth. The latter may be approximated by a theta function in 



momentum space, such as |23 



exp 



/ 5m^ \ ^ 
-4.3 (2.25) X 10 V -^ 



(65) 



derived from Eq. ( p3|) for a = e(/i, r). Substituting Eq.(|6^) into the above, one finds 
that DHPS's upper bound corresponds to having a final sterile neutrino distribution of 
z'^ = 0.04 (0.02). DHPS found an expression similar to Eq. (^) in the limit s <^ 1 for a = e 
[see their Eq. (54)]. They argue that such accumulation of sterile neutrinos is able to destroy 
completely the asymmetry generation since the factor 1 — zf appearing in the term A in the 
growth rate [Eq. ^ now favors stability. 

As explained before, in the regime prior to the onset of asymmetry generation, the 
static approximation and DHPS evolution equations are in complete agreement at leading 
order; this suppression effect is to be distinguished from DHPS's supposedly new but as 
yet negligible "back-reaction". On the other hand, any interested person can numerically 
integrate the static approximation equations at leading order and find no abrupt cut-off 
even when the sterile neutrino distribution increases to the level of ^ 0.04. Therefore 
the severity of DHPS's upper bound has to be imputed to numerical errors.^ 

Physically, the negative contribution from a non-negligible sterile neutrino distribution 
associated with a large sin 26 is compensated for by an increase in the critical momentum. 
In more descriptive words, the asymmetry generation "waits" for the onset of instability 
when this compensation occurs at lower temperatures (see Ref. [ESI for further discussions). 



C. Chemical potentials and final value of neutrino asymmetry 

As noted in the Sec. pi C[ the role of chemical potentials prior to the onset of explosive 
growth is to prevent the asymmetry from being diminished to arbitrarily small values. This 
function is described by the term B in the static approximation [Eq. (pQ])], or, equivalently, 
contained in the variable s{za — Za) (or Fa- in DHPS's notation) in the QKEs [Eq. (|5D|)], 
which DHPS ignored. Another mechanism that can inhibit drastic lepton number destruc- 
tion is the presence of a sterile neutrino asymmetry, created by oscillations from an initial 



The calculation of zf should be done numerically, coupled with the asymmetry evolution equa- 
tion. The analytical expression in Eq. ( |65| ) approximates the production with a theta function 
which is actually a sharp jump around the resonance, valid only before the onset of neutrino 
asymmetry generation (see Ref. [^ for more details). 
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a-neutrino asymmetry. This effect is manifested in a nonzero s{zs — Zg) (or Fs-) in the 
QKEs, and corresponds to the term B' in the static approximation. In the case of a zero 
initial asymmetry, this term is neghgible. On the other hand, the term B due to finite chem- 
ical potentials must never be neglected. The fact that DHPS omit both B and B' causes 
their solution to collapse to values of order 10^^°° prior to the onset of asymmetry growth. 

Furthermore, observe that the condition of lepton number conservation [Eq. (|30D ] and 
the expressions for Pq and Pq in Eq. together lead to the following corollary: 

m - - i^a - za)] f:^{y)y'dy ^ 0. (66) 

Clearly, if one neglects the chemical potential /Iq,, the remaining integral is no longer van- 
ishing, but becomes a fictitious friction force proportional to — (E) L,^^ in Eq. (^) through 
the otherwise absent — term. This point was noted in Ref. |jT9| and speculated to 



be culpable for DHPS's distorted results. However (and we are arguing in DHPS's favour), 
this phantom force does not in fact appear in DHPS's final evolution equation (51), as they 
specifically define J^T{y) {z^ — Za)f°^{ii)y'^dy ^ to be their lepton charge conservation 
condition [their Eq. (30)] which, in the absence of chemical potentials, is actually mathe- 
matically consistent with Eq. (PP]). In effect, DHPS have taken, though fortuitously, the 
appropriate measures to negate the artificial friction force inherent in their formulation. 



VI. CONCLUSIONS 

We review how the numerical solution to the QKEs supplies a lepton asymmetry growth 
curve whose features can be physically understood and reproduced by approximate evolution 
equations. In the higher temperature collision dominated epoch, the static approximation 
mimics the explosive amplification extremely accurately, and allows one to understand why 
a critical temperature exists for 5m? < and cos26'o ~ 1. As the collision rate decreases 
with temperature, the MSW effect eventually takes over as the dominant amplifier, since 
the transitions are adiabatic for a large range of mixing parameters. Also, the separation 
of the neutrino and antineutrino resonance momenta after a significant L has been created 
allows efficient conversion of antineutrinos (for L > 0) into sterile states, while oscillations 
are matter-suppressed for neutrinos (and vice versa for L < 0). The roles of the static 
approximation and the MSW effect are succinctly depicted in Figs, la and 2a. 

The major flaw in the DHPS analysis is their neglect of the pure MSW effect in the post 
explosive growth epoch. Their Eqs. (42) and (51) are approximations to the QKEs valid in 
the collision dominated regime and for a sufficiently small L. Like the static approximation, 
they cannot describe the MSW effect, and therefore underrate the asymmetry growth at 
lower temperatures. However, the DHPS underestimation is much more severe than the 
static case, because their Eq. (51) contains a bogus back-reaction term Bi that becomes 
large and artificially terminates the asymmetry growth under precisely the same conditions 
that imply MSW dominance. This can be understood from the symbolic relation between 
the two approximation schemes: (Y+x) — 1 where the left hand side represents the static 
situation and the right hand side correlates with DHPS. The back reaction term x leads to 
a severe, but fictitious, cut-off when the MSW effect dominates at x ~ 1. 
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We conclude that origin of the 0.2 — 0.37 range for the final steady state value of the 
asymmetry is well understood analytically and physically, in a way which is completely 
consistent with brute force numerical solutions of the Quantum Kinetic Equations. 
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APPENDIX A: TRANSLATION OF DHPS NOTATION 

We give conversion relations between the DHPS nomenclature and the notation in this 
Comment, denoted by the sign with DHPS on the left and our equivalent on the right. 
The effective potential is 

K"// = ±C,rjGFT' + C^^i^ - -K, 



C, = ^^^^ 0.345, = 2v/2C(3) >1 ^ ^ ^^^^ ^^ .^g^s), (Al) 



for a = e (/i, r), and E is the neutrino energy. The quantity 

77 = L^'^) (A2) 

is the so-called effective asymmetry (abbreviated to L in the body of this Comment). 
The neutrino matrix and its expansion in terms of the Pauli matrices are 

r,(r,\ — ( P"-'^ Pas \ ^ fcq [ Pq + Pz Px " iPy \ / a q\ 

^^^^ " \Psa Pss)~ 2 [P. + tPy Po-Pj- ^^"^^ 

For the four matrix elements, DHPS write down the following kinetic equations: 

iPaa = Fo{psa - Pas)/2 - iTo{paa - feq), 
ipss = -Fo{psa - Pas)/'2, 

iPas = WoPas + FoiPss - Paa) / - iPas, 

iPsa = -WoPsa - Fo{pss - Paa)/2 - iTiPsa, (A4) 

with the quantities Wq and Fq given by the expressions 

Fo = —- sin 200 = V,, Wo = —- cos 200 + K// = -V,, (A5) 
and we have used E p. The rates Vq and Fi for a = e (/x, r), and their translations are 
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To = 2T, = ga^-^^GlT'-^ = 1.13 (0.794)G^T^y, 
To = r, Ti = D. (A6) 



However, the coefficient ga/S.15 is smaller than ka [see Eq. (|3^) ] found in Ref. [|T5|. This 
discrepancy simply leads to a different correspondence of the final asymmetry with a given 
choice of parameters (sin^ 29q, 6171^), and is thus irrelevant for the solution of the controversy. 

Are the two sets of kinetic equations (|^) and ( |A^ ) equivalent? Rewriting the former in 
terms of Px,y,z,o by way of Eq. ( |A3|) , one finds that they are indeed concordant, save for 
the replacement Vz —>■ —Vz which is just a matter of convention. Thus we conclude that 
approximation schemes in DHPS and in Refs. /(j|,|^[71|/ begin with identical QKEs. However, 
as discussed in the main text, DHPS use the Fermi-Dirac distribution with zero chemical 
potential as the active neutrino equilibrium distribution function. 

Moving on, DHPS replace the variables pij with ones closely connected to Pq and P: 

Paa = /eq[l + o] ^ = ^(Pq + P.) " 1 = 2^ " 1, 
Pss ^ + S]^S= i(Po - P,) - 1 = 2;, - 1, 

Po.s^f.^[h + ^l]^h=^, 1 = -^. (AT) 

They then take sums and differences of corresponding functions for neutrinos and antineu- 
trinos, i.e., for any generic variable X and its antineutrino counterpart X, they define 

^ (A8) 

The time integration variable is changed to the inverse of the dimensionless temperature 

MeV , 
X = -j^, (A9) 

which generates a factor Tix. Therefore DHPS introduce the following new quantities: 

= H = £- 

TCx TCx^ 

Hx 2pTi.x Hx ' 

Hx 2pHx ° ~^ Hx Hx ^ ^ 

They also find it useful to separate W in two different contributions, W = U ± VZ, with 

VZ=^-^I^^-^a, (All) 

Hx 2pHx 
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Z =10^' Lo- r%'f>-]- (A12) 



where 



\ Jo Att^ 

The quantity Z, as defined above, is clearly connected to the lepton asymmetry r], although 
their exact relationship is not specified by DHPS. This we provide here for completeness: 

With these definitions, DHPS rewrite the QKEs in the new variables s±, a±, h±, l±: 



4 


= Fl±, 






«± 


= -Fl± 


- 27+a± - 


27^a^, 




= Ul± - 




-h± - 7-^T' 



4 = f («± - s±) - Uh± + VZh^ - 7+/± - 7_/:p, (A14) 

in which terms containing 7_ = (Fi — Ti)/2 are later discarded. This is also the case in Eq. 
(I). Thus we stress again that the DHPS kinetic equations are equivalent to Eq. (H). 

Now DHPS begin to extract their evolution equations by seeking approximate solutions 
to Eq. ( |A14{ ). More new notations and symbols are introduced during the derivation, in 
particular, the integration variable x is replaced with r by way of the function U: 

f/ = cos 2^07^ f^-lV (A15) 



2pHx \t 

such that U = at T = 1. Then it follows that 



where ^ ~ 6.5 (12.8) x lO^cos 2^0^, for a = e (/i, r). 

At this point DHPS divide all the equations by the quantity M, where 

M=^^1.12xl0^^\^, (A17) 

and make the replacements ]g(7, F, U, VZ) — > (7, F, f/, VZ), leading to 





2p r _ ^ ^ 


7 = 




F = 


— tan2^^ ~ —s, 






U = 






r"^ c 



7^0 

\ res / 

a 



- 1, 



l^Z = - ~ a. (A18) 
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These new definitions combine to produce the change 



d 1 d 
dx Q dr^ 

in the left hand side of the kinetic equations, where 

-1 



(A19) 



COS 26*0. 



(A20) 



More new definitions follow from here, whose translations are hsted below (with c ~ 1): 
D^ZV = a, 

a^^^^ + U^^d^ + {c-bf, 

/ D2 



^(6-1) 



1 - 



(1 -6)2 

,2 



1 + 



d^ + {l- 6)2 



1 



t = - 



T 



- ~ Vb. 



(A21) 



The quantity r is inversely proportional to y. DHPS thus define the new variable 



'rpO I 

T 



(A22) 



Moreover they introduce a new function 60 that expresses the difference between sum dis- 
tribution of active neutrinos and sum distribution of sterile neutrinos: 



60 ~ -s+ = 1-2;+, 

such that their final sterile neutrino production equation [Eq. (42) in DHPS] is 

dbo _ QF^^, 

while the asymmetry evolution equation [their Eq. (51)] is given by 



(A23) 



(A24) 



t^(f - l)/e,(tg)6o(l/t) 



(72^2 



F2(a2 + 

4(72^2 



(A25) 



With some trivial algebra one can show that 



^2^2 = [72 + ([/ + £,)2][^2 ^ _ ^)2]^ 
0-2 + ^2 ^ ^2 ^ ^2 ^ 



(A26) 
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Using the conversion relations, the DHPS sterile neutrino production equation translates to 



9(T/MeV) mr A(0) '^^ 

while their neutrino asymmetry evolution equation is equivalent to 



dL 



A 



1 - 



g2 d^ + (c- hf + 
4 A(0) 



(A27) 



(A28) 



Lrf(T/MeV) 7o A(0) 
in our language. 

Note Added 



Shortly after completion of this paper, R. Buras and D. V. Semikoz informed us of their work 
( |hep-ph/0008263|) , where they also concluded that the main result of DHPS was incorrect. 
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FIGURES 



FIG. 1. Evolution of the effective total lepton number, L^'^\ for i/j- Vg oscillations for the 
parameter choice, Sm? = — lOeV^, sin^2^o = 10~^. In Fig. la, the solid curve is the numerical 
solution of the QKEs while the dashed curve is the static approximation for rlrcs = (^int/^m)|rcs < 1 
while for r|rcs > 1 pure adiabatic MSW evolution is employed. The dashed dotted line is the static 
approximation (over all T). Fig. lb contains the corresponding evolution of the MSW resonance 
momentum/T (i.e. j/res) for the neutrinos (solid line) and anti-neutrinos (dashed line) for this 
example. 

FIG. 2. Same as Fig. 1 except for the parameter choice 5m^ = — lOOeV'^ and sin^ 2^o = 10^^. 
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